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Circulating Muscle-specific miRNAs in Duchenne 
Muscular Dystrophy Patients 

Xihua Li 1 , Yuying LP, Lei Zhao 1 , Duo Zhang 2 , Xuan Yao 2 , Huihui Zhang 2 , Yu-cheng Wang 34 , Xin-yi Wang 4 , Hongfeng Xia 2 3 , Jun Yan 5 and 
Hao Ying 23 

Noninvasive biomarkers with diagnostic value and prognostic applications have long been desired to replace muscle biopsy for 
Duchenne muscular dystrophy (DMD) patients. Growing evidence indicates that circulating microRNAs are biomarkers to assess 
pathophysiological status. Here, we show that the serum levels of six muscle-specific miRNAs (miR-1/206/133/499/208a/208b, 
also known as myomiRs) were all elevated in DMD patients (P < 0.01). The receiver operating characteristic curves of circulating 
miR-206, miR-499, miR-208b, and miR-133 levels reflected strong separation between Becker's muscular dystrophy (BMD) and 
DMD patients (P < 0.05). miR-206, miR-499, and miR-208b levels were positively correlated with both age and type lie muscle 
fiber content in DMD patients (2-6 years), indicating that they might represent the stage of disease as well as the process of 
regeneration. miR-499 and miR-208b levels were correlated with slow and fast fiber content and might reflect the ratio of slow to 
fast fibers in DMD patient (>6 years). Fibroblast growth factor, transforming growth factor-p\ and tumor necrosis factor-a could 
affect the secretion of myomiRs, suggesting that circulating myomiRs might reflect the effects of cytokines and growth factors 
on degenerating and regenerating muscles. Collectively, our data indicated that circulating myomiRs could serve as promising 
biomarkers for DMD diagnosis and disease progression. 
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Introduction 

Duchenne muscular dystrophy (DMD) is a lethal X-linked 
recessive neuromuscular disease affecting 1 in 3,500 new- 
born males. DMD is caused by mutations in dystrophin, a 
cytoskeletal protein that connects the cytoskeleton of a mus- 
cle fiber to the surrounding extracellular matrix through the 
cell membrane. Mutations leading to complete absence of 
dystrophin protein cause myofiber necrosis, accompanied by 
a complex series of events including phagocytosis, infiltra- 
tion of inflammatory cells and decrease of muscle fibers with 
subsequent fibrosis and fat replacement leading to subse- 
quent loss of ambulation and premature death by respiratory 
or heart failure. 12 Loss of dystrophin renders myofibers sus- 
ceptible to exercise induced injury, which triggers continuous 
degeneration and regeneration. 3 In contrast, mutations lead- 
ing to the production of less defective yet still partially func- 
tional dystrophin protein result in a much milder dystrophic 
phenotype in affected patients. These patients were known 
as Becker's muscular dystrophy (BMD) patients. 4 

As for patients with suspected DMD, the preliminary 
diagnostic criterion is the serum creatine kinase (CK) level, 
which is markedly increased in early stages of the disease. 
However, it is not always reliable since CK levels increase in 



normal individuals due to exercise or muscle trauma. More- 
over, it is not suitable as a screening test for those who have 
already been using wheelchairs, because CK levels may fall 
when disease progresses, probably due to the loss of muscle 
tissue. 5,6 In addition, the CK level in BMD patients is also 
highly elevated, which makes CK level a poor biomarker to 
distinguish DMD from BMD. The precise diagnosis for DMD 
should contain a combination of genetic analysis, muscle 
biopsy, and clinical observation of muscle strength and 
function. Nonetheless, muscle biopsy is invasive and pro- 
vides only local information, which may not reflect the state 
of all muscles, while muscle strength assessment is subject 
to extensive intra- and inter-patient variability. 7 Therefore, 
more reliable biomarkers with diagnostic value and prog- 
nostic applications have long been desired. Although there 
is no cure for DMD, several strategies including gene ther- 
apy, cell therapy, and pharmacotherapy have been recently 
tested in animal models of DMD 89 and some of them have 
already entered clinical trials. Considering the various clinical 
treatments and trials currently ongoing for DMD, there is an 
urgent need to develop better noninvasive biomarkers that 
could be used to monitor the progression of DMD pathology, 
guide therapeutic decisions, and evaluate the efficacy of 
potential therapies. 
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MicroRNAs (miRNAs) are important regulators for numer- 
ous physiological and pathological processes. 10-12 Growing 
evidence revealed that miRNAs also exist in various body 
fluids such as serum, plasma, saliva, urine, and milk. 13-15 
Altered profiles of circulating miRNAs have been shown 
to associate with various diseases, including tissue injury, 
cancers, and diabetes. Consequently, circulating miRNAs 
become promising biomarkers to assess the pathophysiolog- 
ical status and monitor the effectiveness of clinical treatment. 
It has been shown that serum levels of several myomiRs 
(miR-1, miR-133a, and miR-206) are increased not only in 
animal models, such as mdx mice and CXMDJ dogs, but also 
in DMD patients. 16-18 In the human study, the authors showed 
that miR-1, miR-1 33, and miR-206 were able to discriminate 
DMD patients from healthy controls and BMD patients. Most 
importantly, the increased levels of these miRNAs but not the 
CK was correlates with the severity of muscle damage show 
a very clear inverse correlation with North Star Ambulatory 
Assessment (NSAA) scores. However, only 26 DMD and 5 
BMD patients were recruited for serum miRNA level study 
and only data from 10 DMD patients between 3 and 6 years 
of age were used for evaluating the correlation between 
serum miRNA levels and clinical assessments. Therefore, 
more comprehensive investigation is required to confirm 
these findings in cohort with more patients across broader 
age groups. 

The major goal of this study was to identify potential circu- 
lating miRNA biomarker for DMD diagnosis and prognosis. 
MyomiRs are suggested to have a role either in myoblast pro- 
liferation, differentiation, and muscle regeneration, which are 
key events in the pathogenesis of DMD. Besides three classical 
myomiRs (miR-1 , miR-1 33, and miR-206), we also investigated 
the serum levels of another two myomiRs, miR-499 and miR- 
208, which have been shown also play critical roles in muscle 
biology and exercise adaptation. 19 miR-499 is expressed mainly 
in the slow skeletal muscle and cardiac muscle; miR-208a is 
expressed mainly in the cardiac muscle and miR-208b in the 
slow skeletal muscle. 2021 We hypothesize that these myomiR 
are indicative for the disease pathophysiology and have poten- 
tial as prognostic biomarkers. We examined serum levels of 
myomiRs in mdx mice, DMD and BMD patients and compared 
the serum levels of miRNAs with clinical assessment including 
age, CK value, and muscle fiber composition. Based on our 
results, we proposed that specific serum circulating miRNAs 
are not only novel biomarkers for DMD diagnosis but also may 
be useful to monitor the pathological progression. 

Results 

Serum levels of myomiRs in dystrophic subjects 

To test the value of serum miRNAs as potential biomark- 
ers for DMD diagnosis, we quantified the levels of myomiRs 
(miR-1, miR-1 33, miR-206, miR-208a, miR-208b, and miR- 
499) in the serum of dystrophin-deficient mice (mdx mice) by 
real-time PCR (Supplementary Figure S1). Consistent with 
an earlier report, 16 the levels of miR-1 , miR-1 33, and miR-206 
were all significantly elevated. We also observed a significant 
increase in the levels of miR-499 and miR-208b. In contrast, 
the serum level of miR-208a did not change in mdx mice. 

Subsequently, we determined the levels of these myomiRs 
in serum samples from healthy (n = 23), Becker (n = 15), 
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Figure 1 Elevation of myomiRs in dystrophic children serum. 

(a-f) Serum levels of miRNAs (miR-1 , miR-206, miR-1 33, miR-499, 
miR-208a, and miR-208b) in dystrophic children and normal controls 
(normal, n = 23; BMD, n = 1 5; DMD, n = 52) were determined by real- 
time PCR. ANOVA was used for statistical analysis. SEs are shown 
as error bars. *P< 0.05, **P< 0.01, "*P< 0.001. BMD, Becker's 
muscular dystrophy; DMD, Duchenne muscular dystrophy; NS, not 
significant; SE, standard error. 



and Duchenne (n = 52) children, aged from 1 to 14 years 
old (Figure 1a-f). These patients were not on glucocorticoid 
treatment and were still ambulant except one patient (Sup- 
plementary Table S1). As expected, the levels of miR-1, 
miR-1 33, miR-206, miR-208b, and miR-499 were signifi- 
cantly elevated in DMD patients compared to healthy control 
subjects. To be noted, the level of miR-208a was also sig- 
nificantly increased in DMD subjects. In addition, the levels 
of these six myomiRs were also elevated in BMD patients, 
although to a less extent and not statistically significant 
except for miR-499. Importantly, the levels of miR-206 and 
miR-499 were significantly higher in DMD than those in BMD 
patients. The levels of miR-1 miR-1 33, miR-208a, and miR- 
208b showed a tendency to be higher in DMD than those in 
BMD; however, the differences did not reach statistical sig- 
nificance. We also compared the CK value in these BMD and 
DMD patients. We observed that the CK values were slightly 
higher in DMD patients than those in BMD but not signifi- 
cantly (Supplementary Figure S2). 

The statistical significance of these measurements and diag- 
nostic power of serum miRNAs were evaluated by receiver 
operating characteristics (ROC) analysis (Figure 2a-f and 
Table 1). Our results demonstrated that miR-206, miR-499, 
and miR-208b displayed better sensitivity and specificity to dis- 
criminate DMD from normal subjects (AUC: 0.9156, 0.99, and 
0.9323, respectively, P < 0.0001; Figure 2b,d,f and Table 1). 



Molecular Therapy — Nucleic Acids 



Circulating myomiRs in DMD Patients 

Lietal. 



3 miR-1 b miR-206 




100% - specificity% 100% - specificity% 

Figure 2 ROC analysis of serum levels of myomiRs. (a-f) ROC 
curves were plotted to evaluate the diagnostic power of myomiRs 
(miR-1, miR-206, miR-1 33, miR-499, miR-208a, and miR-208b). 
Red line, DMD versus normal; blue line, BMD versus normal; green 
line, BMD versus DMD. The area under the ROC curve (AUC) was 
calculated for the measurement of discrimination accuracy (normal, 
n = 23; BMD, n = 15; DMD, n = 52). BMD, Becker's muscular 
dystrophy; DMD, Duchenne muscular dystrophy; ROC, receiver- 
operator characteristic. 

Also miR-1 , miR-1 33, and miR-208a displayed a very good sta- 
tistical scores (AUC: 0.8227, 0.8119, and 0.8127 respectively, 
P < 0.0001; Figure 2a,c,e and Table 1). Similarly, miR-206, 
miR-499, and miR-208b had higher sensitivity and specificity 
to discriminate BMD from controls (AUC: 0.8725, 0.9913, and 
0.913, respectively, P< 0.0001), while miR-1, miR-1 33, and 
miR-208a were less potent (AUC: 0.8406, 0.7333, and 0.8261 , 
respectively, P= 0.0005, P= 0.0163, and P= 0.0008, respec- 
tively; Figure 2a-f and Table 1). Most importantly, ROC curve 
analysis displayed the diagnostic power of serum miR-208b, 
miR-206, miR-499, and miR-1 33 in discriminating DMD from 
BMD patients (AUC: 0.71 1 5, 0.7090, 0.6987, and 0.6756, listed 
from high to low, P = 0.0131, P = 0.0143, P = 0.0198, and 
P= 0.0394, respectively; Figure 2f,b,d,c and Table 1). We also 
performed ROC analysis to evaluate the diagnostic power of CK 
value in discriminating DMD from BMD. As shown in Table 1, 
CK did not have diagnostic power in discriminating DMD from 
BMD patients (AUC: 0.6536, P = 0.0723). These results indi- 
cated that circulating miRNA signatures have diagnostic value 
in patients with DMD. 

The correlation between serum level of myomiR and age 

The absence of dystrophin results in progressive degenera- 
tion of skeletal and cardiac muscle. We speculated that the 
circulating miRNA levels might change during aging and 
might correlate with the extent of muscle degeneration or 



Table 1 Receiver-operator characteristic (ROC) curve analysis of miRNAs 
and CK for diagnosis of dystrophic children 





Normal vs DMD 


Normal vs BMD 


BMD vs DMD 


AUC 


P value 


AUC 


P value 


AUC 


P value 


miR-1 


0.8227 


<0.0001 


0.8406 


0.0005 


0.5769 


0.3668 


miR-206 


0.9156 


<0.0001 


0.8725 


0.0001 


0.7090 


0.0143 


miR-133 


0.8119 


<0.0001 


0.7333 


0.0163 


0.6756 


0.0394 


miR-499 


0.9900 


<0.0001 


0.9913 


<0.0001 


0.6987 


0.0198 


miR-208a 


0.8127 


<0.0001 


0.8261 


0.0008 


0.5603 


0.4796 


miR-208b 


0.9323 


<0.0001 


0.9130 


<0.0001 


0.7115 


0.0131 


CK 










0.6536 


0.0723 



BMD, Becker's muscular dystrophy; CK, serum creatine kinase; 
DMD, Duchenne muscular dystrophy. 
Normal, n = 23; BMD, n = 1 5; DMD, n = 52. 



the stage of the disease. To test our hypothesis, we first per- 
formed a case-by-case correlation analysis between age and 
serum miRNA/CK levels of dystrophic children of all ages. 
As expected, we did not observe any statistical significance 
of correlations (Supplementary Table S2). Subsequently, 
we analyzed the correlation of serum miRNA or CK levels 
with age in different age groups (<2, 2-6, and >6 years old; 
Supplementary Table S2). Interestingly, serum miR-206, 
miR-499, and miR-208b levels were positively correlated 
with age in DMD patients aged 2-6 years (Figure 3a-c and 
Supplementary Table S2). No correlation between CK value 
and age was observed in DMD or BMD patients in any age 
groups (Supplementary Table S2). These results indicated 
that these age-associated serum miRNA signatures might 
represent the extent of muscle degeneration or the stage of 
the disease. 

The correlation between serum miRNA level and muscle 
fiber type composition 

Since myomiRs play regulatory roles in muscle growth, dif- 
ferentiation, fiber determination, and metabolism, we hypoth- 
esized that the serum levels of these miRNA might reflect the 
muscle structure and function and correlate with the patho- 
genesis of DMD or BMD. In order to test our hypothesis, we 
first determined the fiber type composition in the muscle- 
biopsy specimens from these DMD and BMD subjects. We 
were not able to observe any correlation of fiber type com- 
position with age in dystrophic subjects of all ages (Supple- 
mentary Table S3). In agreement with the concept that fast 
muscle fibers are preferentially affected, 22 we found that the 
percentage of fast-twitch (type Ma and type lib) muscle fibers 
decreased in DMD patients aged 2-6 years (Figure 3d and 
Supplementary Table S3). In contrast, in the same age 
group (2-6 years), the percentage of type lie muscle fibers 
increased gradually (Figure 3e and Supplementary Table 
S3), while in patients older than 6 years the percentage of 
type lie muscle fiber decreased (Figure 3f). Because many 
of type lie muscle fibers have characteristics of regenerat- 
ing fibers and reflect an active regenerating process, 23-26 our 
result suggested that the regeneration becomes increasingly 
active in patients between the ages of 2 and 6 years due to 
continuous degeneration, while in older patients (>6 years) 
after degeneration and regeneration cycle depletes the sat- 
ellite cell pool, these patients fail to regenerate and recover 
muscle mass. We did not observe any correlation between 
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type lie fibers and age in all age groups of BMD, in turn, sug- 
gesting that muscle regeneration in BMD patients remains 
active (Supplementary Table S3). 

Then, we evaluated the correlation between serum miRNA 
levels and muscle fiber composition. As shown in Table 2, 
serum level of miR-499 was negatively correlated with fast- 
twitch fiber composition, while serum level of miR-133 was 
positively correlated with fast-twitch fiber composition and 



negatively correlated with slow-twitch fiber composition 
(Table 2). Since both serum miRNA levels and muscle fiber 
composition might be correlated with age at different stages 
of disease, in order to obtain exact correlations, we ana- 
lyzed the correlation between them in different age groups. 
We found that the serum levels of miR-206, miR-499, and 
miR-208b were positively correlated with the percentage of 
type lie muscle fibers in DMD patients between 2 and 6 years 
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Figure 3 Correlation between serum miRNA level or muscle fiber type composition and age. (a-c) Case-by- case comparison of miR- 
206 (a), miR-499 (b), and miR-208b (c) levels with age in 20 DMD children between 2 and 6 years of age. Regression lines are displayed. 
(d,e) The percentage of fast-twitch (type lla and type lib) muscle fibers (d) shows an inverse correlation, while the percentage of type He (e) 
shows a positive correlation with age in DMD children from 2 to 6 years old (n = 1 5). (f) In DMD patients older than 6 years, the percentage of 
type lie shows a negative correlation with age. DMD, Duchenne muscular dystrophy; r, correlation coefficient. 



Table 2 Correlation between muscle fiber type composition and serum miRNAs/CK level of DMD children 



Muscle type 




<2 years (n = 6) 


2-6 years (n = 12) 




>6 years (n = 1 1 ) 


0-14 years (n 


= 29) 


(Type lla + llb)% 


r 


P value 


r 


P value 


R 


P value 


r 


P value 


miR-1 


0.0320 


0.9520 


0.1463 


0.6501 


-0.3084 


0.3561 


0.0829 


0.6688 


miR-206 


-0.2437 


0.6416 


-0.4142 


0.1807 


-0.5872 


0.0575 


-0.1736 


0.3678 


miR-133 


0.7682 


0.0744 


0.2332 


0.4657 


-0.5159 


0.1043 


0.4212 


0.0229 


miR-499 


-0.0974 


0.8544 


-0.5610 


0.0578 


-0.6931 


0.0180 


-0.3771 


0.0437 


miR-208a 


0.5516 


0.2565 


0.2731 


0.3904 


-0.2957 


0.3773 


0.1489 


0.4406 


miR-208b 


0.7433 


0.0904 


-0.1409 


0.6624 


-0.6153 


0.0439 


0.3609 


0.0544 


CK 


0.0805 


0.8795 


-0.3846 


0.2170 


-0.0801 


0.8149 


-0.1991 


0.3006 


Type 1% 


r 


P value 


r 


P value 


R 


P value 


r 


P value 


miR-1 


-0.3730 


0.4665 


-0.5457 


0.0665 


0.4052 


0.2163 


-0.3373 


0.0736 


miR-206 


0.1673 


0.7514 


-0.1964 


0.5407 


0.6807 


0.021 1 


-0.0470 


0.8086 


miR-133 


-0.6527 


0.1600 


-0.4771 


0.1168 


0.6136 


0.0447 


-0.4048 


0.0294 


miR-499 


0.0865 


0.8706 


0.0126 


0.9690 


0.7047 


0.0155 


0.1666 


0.3877 


miR-208a 


-0.3631 


0.4793 


-0.3064 


0.3327 


0.3775 


0.2524 


-0.1304 


0.5001 


miR-208b 


-0.5348 


0.2742 


-0.3460 


0.2706 


0.6856 


0.0199 


-0.3480 


0.0643 


CK 


-0.1725 


0.7438 


0.2768 


0.3837 


0.1745 


0.6078 


0.2155 


0.2614 


Type llc% 


r 


P value 


r 


P value 


R 


P value 


r 


P value 


miR-1 


0.2178 


0.6784 


0.5474 


0.0654 


-0.1602 


0.6379 


0.2957 


0.1194 


miR-206 


0.0237 


0.9644 


0.7201 


0.0083 


-0.1073 


0.7534 


0.2382 


0.2133 


miR-133 


-0.4398 


0.3829 


0.3606 


0.2495 


-0.1271 


0.7095 


-0.0852 


0.6605 


miR-499 


-0.1305 


0.8054 


0.6122 


0.0344 


0.0874 


0.7984 


0.2639 


0.1666 


miR-208a 


-0.4243 


0.4018 


0.0934 


0.7729 


-0.1293 


0.7047 


-0.0371 


0.8486 


miR-208b 


-0.5803 


0.2273 


0.6088 


0.0356 


-0.0518 


0.8797 


-0.0954 


0.6224 


CK 


-0.0195 


0.9707 


0.0701 


0.8286 


-0.1916 


0.5726 


-0.0155 


0.9366 



CK, serum creatine kinase; DMD, Duchenne muscular dystrophy. 
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Figure 4 Correlation between serum miRNA level and type lie or type I muscle fiber composition, (a-c) Serum levels of miR-206 (a), miR-499 
(b) and miR-208b (c) show positive correlation with type He muscle fiber composition in DMD children between 2 and 6 years old (24< age 
<72; n = 12). (d-f) Serum levels of miR-206 (d), miR-499 (e) and miR-208b (f) show positive correlation with type I muscle fiber composition 
in DMD children older than 6 years old (age > 72; n- 1 1). DMD, Duchenne muscular dystrophy; r, correlation coefficient. 



old (Figure 4a-c and Table 2), suggesting that the serum 
miR-206, miR-499, and miR-208b might serve as real-time 
biomarkers to monitor the activity of regeneration in DMD 
patients between 2 and 6 years old. In addition, we observed 
that serum level of miR-499 and miR-208b were positively 
correlated with the percentage of slow-twitch (type I) mus- 
cle fibers and negatively correlated with fast-twitch (type Ma 
and lib) muscle fiber composition in DMD patient older than 
6 years (Table 2). Similar correlations were observed for 
miR-206 although one of them tended to be significant (P = 
0.0575; Table 2). This observation suggested that the serum 
levels of miR-206, miR-499, and miR-208b might reflect the 
ratio of slow to fast muscle fibers in DMD patients older than 6 
years (Figure 4d-f). Interestingly, these three miRNAs were 
mainly expressed in solues muscle (slow-twitch muscle) of 
mice (Supplementary Figure S3). In contrast to these three 
miRNAs, miR-1, and miR-133 were enriched in tibialis ante- 
rior muscle (TA, fast-twitch muscle) of mice (data not shown). 
These observed correlations between serum miRNA levels 
and muscle fiber composition might be due to the differential 
expression pattern of miRNA in muscles comprising of differ- 
ent fiber types. As we expected, no correlation between CK 
value and muscle fiber composition was observed in all age 
groups (Table 2). 

The regulation of miRNA secretion by cytokines and 
growth factors 

The exploration of circulating miRNAs brings about an 
important advance for biomarker discovery. Circulating 
miRNAs have been profiled in many diseases. However, 
the source of these circulating miRNAs and the underlying 
mechanism are not fully elucidated. We first tested whether 
the serum from DMD patient could affect the secretion of 
myomiRs. Interestingly, we found that DMD serum increased 
the miR-1, miR-133, miR-206, miR-499, and miR-208b lev- 
els in the medium of differentiated C2C12 cells (Figure 5a). 



We also directly measure the miRNA levels with the same 
amount of the DMD serum in the same amount of medium 
and detected very low levels of these miRNAs, suggesting 
that the human miRNA contamination from added DMD 
serum is limited (Figure 5a). 

It is known that cytokines and growth factors are enriched 
in regenerating muscles and play important roles in inflam- 
mation and regeneration. 2728 To figure out which factors or 
molecules in the DMD serum might contribute to miRNA 
secretion, we investigated the effect of tumor necrosis 
factor-a (TNF-a), fibroblast growth factor (FGF), and trans- 
forming growth factor-p 1 (TGF-P). As shown in Figure 5b, 
these factors have differential impact on the level of myomiRs 
in the culture medium of C2C12 myotubes. Moreover, consis- 
tent with the finding that FGF could increase the level of these 
myomiRs in the culture medium, inhibition of FGF signaling 
by PD1 73074, a specific inhibitor for fibroblast growth factor 
receptor 1 (FGFR1), resulted in a reduction of myomiR levels 
in the culture medium (Supplementary Figure S4). Further- 
more, the effect of FGF could be blocked by pre-incubating 
C2C12 myotubes with PD1 73074 (Figure 5c). In addition, 
TGF-p affected the level of miR-1, miR-133, miR-206, and 
miR-208b but had little effect on miR-499 level, while TNF-a 
increased the level of miR-133, miR-206, and miR-208b, but 
not miR-1 and miR-499 (Figure 5b). To rule out the possibil- 
ity that treatment of TNF-a, FGF, and TGF-p could disturb the 
membrane integrity, leading to the leak of some miRNA out 
of the cells, we checked the effect of these growth factors and 
cytokines on the cell viability. As shown in Supplementary 
Figure S5a-d, the treatment did not significantly affect the 
cell viability of C2C12 myotubes. 

Based on these findings, we proposed that the differ- 
ent circulating miRNA profiles might also reflect the action 
of different cytokines and growth factors on degenerating 
and regenerating muscles. On the other hand, the differen- 
tial effects of cytokine and growth factors on the secretion 
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Figure 5 The effect of DMD serum, cytokine, and growth factors on the level of myomiRs in culture medium of C2C12 myotubes.The levels 
of myomiR (miR-1 , miR-1 33, miR-206, miR499, and miR-208b) in the culture medium of C2C1 2 myotubes were examined by using real-time 
PCR. (a) C2C12 myotubes were treated with 20 ul serum from DMD patients or control subjects. The culture medium from these two groups 
were collected for measurement and designated as medium-DMD serum and medium-Control serum, respectively, as indicated. Since 
the primers could not distinguish human and mouse miRNAs, a third group was set to measure the miRNA levels in 20 ul DMD serum and 
designated as DMD serum as indicated. Student t test was used to analyze the data from DMD serum treatment and control serum treatment, 
(b) C2C12 myotubes were incubated with FGF, TNF-a, orTGF-fS. (c) Differentiated C2C12 myocytes were preincubated with FGFR1 inhibitor 
(PD1 73074) before FGF was added to the culture medium. Error bars represent the SEs. Student's f test was used. *P< 0.05, **P< 0.01, ***P 
< 0.001 versus control group. DMD, Duchenne muscular dystrophy; SE, standard error. 



of different miRNAs suggested that the secretion might be 
regulated by separated mechanisms, which requires further 
investigation. 

Discussion 

Comprehensive approaches including gene therapy, cell 
therapy, and pharmacotherapy have been proposed for an 
ultimate cure for DMD. Considering the various clinical treat- 
ments and trials currently ongoing for DMD, the identification 
of serum biomarkers can be valuable to diagnose, prognosis 
of diseases, and guide therapeutic decisions. Although CK 
is commonly used as a biomarker of muscular diseases to 
evaluate the level of muscle damage and necrosis, 29 but it is 
not always reliable since CK level varies considerably upon 
several independent stress conditions, such as exercise. 
Growing evidence has indicated that serum miRNAs are 
promising biomarkers for diseases, like tissue injury, cancers, 
and cardiovascular disease. 153031 The complexity of the regu- 
lation and function of miRNA provides us more opportunities 
to identify disease-associated biomarkers. Here our study 
suggested that specific circulating miRNAs in serum could 
serve as potential biomarker for DMD diagnosis. We focused 
on myomiRs since these functional miRNAs have been stud- 
ied extensively in muscles. The elevation of serum levels of 
myomiRs (miR-1, miR-1 33, miR-206, miR-208a, miR-208b, 
and miR-499) were observed in DMD patients. Based on our 
result, among these miRNAs, miR-499 is the best biomarker 
to discriminate DMD from healthy control subjects (AUC = 
0.99). Although miR-208b is not the best biomarker to dis- 
criminate DMD from healthy controls, it displayed the great- 
est power to distinguish DMD and BMD (AUC = 0.71 1 5). Our 



study here suggested that serum miRNA profiles might be 
considered as biomarkers for DMD diagnosis. 

An early report suggested that the levels of serum miRNAs 
might be associated with pathogenesis of DMD. 16 Different 
fiber compositions contribute to muscle functions and features 
such as contraction and metabolism. During degeneration 
and regeneration, muscle undergoes composition changes 
in muscle fiber types. We hypothesized that varied fiber type 
compositions might reflect the structure and function of mus- 
cle and might be used as a signature of disease progression. 
In this study, we compared the muscle fiber composition with 
age and explored the correlation between serum miRNA lev- 
els and muscle fiber composition. We uncovered that serum 
levels of miR-206, miR-499, and miR-208b are positively cor- 
related with age and type lie muscle fibers from 2 to 6 years 
in DMD patients (Figures 3 and 4). Therefore, we suggested 
that these miRNAs might be useful in monitoring the activ- 
ity of regeneration in DMD patients between 2 and 6 years 
old. Moreover, these biomarkers are not replaceable since 
muscle biopsy is invasive, after first diagnosis. 

We also showed that miR-206, miR-499, and miR-208b 
are positively correlated with slow-twitch muscle fibers in 
DMD patients more than 6 years old. It is possible that the 
serum miRNA levels represent the residual muscle mass. 
Hence these results allowed us to propose that the serum 
miRNAs could be valuable prognostic biomarkers to monitor 
DMD diseases progression. This result is consistent with the 
concept that miR-206, miR-208b, and miR-499 are enriched 
in slow-twitch muscle fiber and fast-twitch muscle fibers are 
the first to degenerate. In other words, the less affected slow 
muscle fibers might contribute to the higher level of serum 
miRNA at late stage of DMD. 
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To elucidate the mechanisms underlying circulating 
miRNA profile variation, we incubated the cultured mus- 
cle cells with DMD serum to mimic the microenvironment 
in vivo and found that DMD serum is capable of increasing 
the level of myomiR in the culture medium. This finding sug- 
gested that the myomiRs in serum might be secreted from 
muscle tissues. Also, DMD serum might contain regulatory 
molecules that could regulate miRNA secretion. A previous 
report showed that the expression of cytokines including 
TNF-a was not only upregulated in DMD muscles but also 
increased with age. 32 It is intriguing that we found TNF-a 
as well as other growth factors (FGF and TGF-(3) could 
increase the levels of muscle specific miRNAs in C2C12 
cells under differentiation culture condition. We speculated 
that the increase of these miRNAs in the serum is secreted 
from differentiating muscle cells through the exact mecha- 
nisms remain to be explored. It has been shown that ini- 
tial myofiber damage is exacerbated by the endogenous 
inflammatory response. 33-35 We, therefore proposed that the 
serum levels of these miRNAs may partly reflect the degen- 
eration-induced inflammatory response in patients between 
2 and 6 years old. This would be consistent with the fact that 
there is no significant change in serum revels of some abun- 
dant miRNAs or housekeeping miRNAs, such as miR-16, in 
DMD patients. Although it is likely that the serum myomiRs 
could also be directly leaked from damaged muscle cells, 
we propose that serum miRNA might be mainly derived 
from selective secretion. Further work is required to validate 
our hypothesis in the future. 

A recent study from Dr. Muntoni's group suggested that 
miR-1, miR-133a, and miR-133b showed a similar age 
pattern: higher and gradually decreasing level in patients 
between 4 to 10 years old and considerably lower levels 
of miRNAs in older patients (from 13 to 17 years of age). 18 
Based on the age profile of miR-1 , miR-1 33a, and miR-1 33b, 
the authors speculated that the three dystromirs can be uti- 
lized as biomarkers for the remaining muscle mass in DMD. 
Our results are in general in agreement with this study that 
the increases in miRNA is only significant in young DMD 
patients (<10 years old), whereas the levels of the same 
miRNA start to drop in the older patients. Since the DMD 
patients recruited in our study are much younger than those 
in their study and only one patient had lost ambulation, we 
have more chance to detect possible correlation between 
circulating miRNAs with disease progression in early stage. 
Actually, we found the serum levels of miR-206, miR-499, 
and miR-208b were positively correlated with age in young 
patient (from 2 to 6 years of age). The more elaborated 
parts of our study are the observation that serum levels of 
miRNA is related to muscle fiber type and growth factors 
might influence the secretion/leakage of miRNA from myo- 
cytes. Together, the data from these two studies provided us 
more information about the change of the amount of miR- 
NAs over time. 

Taken together, we identified serum miRNA signatures 
for DMD prognosis. The signatures are associated with 
fiber type composition and disease progression. Thus 
analysis of the pattern of myomiRs could be useful for 
assessing disease severity and effectiveness of therapeu- 
tic interventions. 



Materials and methods 

Mouse and human studies. Animals were maintained and 
experiments were performed according to the protocol 
approved by the Animal Care and Use Committees of Insti- 
tute for Nutritional Sciences (permit number: 2011-AN-14). 
Mice were housed in a temperature controlled room (22 °C) 
with a 12-hour light/dark cycle and had free access to food 
and water. Mdx mice were from Dr.Xiang Gao (Nanjing Uni- 
versity, Nanjing, China). Eight weeks old male mice were 
used in this study. Mouse muscle pads were harvested and 
frozen in liquid nitrogen and stored at -80 °C. Blood samples 
were centrifuged to prepare serum and stored at -80 °C. 

Human study was reviewed and approved by Ethics Com- 
mittee of Children's Hospital of Fudan University. All human 
participants gave written informed consent. The serum sam- 
ples were collected at Children's Hospital of Fudan University. 
Patient inclusion criteria were: genetically and muscle biopsy 
proven DMD or BMD diagnosis, patient still ambulant with- 
out any help (only for children younger than 6-years-old), no 
severe or moderate learning difficulties or behavioral prob- 
lems. The bicep (upper arm muscle) which is about 5mm in 
diameter and 1 cm in length was selected for muscle biopsy 
analysis and myosin ATPase Staining. 

Reagents, cell culture and differentiation. Mouse C2C12 
myoblasts (gift from Dr. Jia Li, SIBS, CAS, China) were 
maintained in a humidified incubator at 37 °C and 5% C0 2 
in Dulbecco's modified Eagle medium containing 10% fetal 
bovine serum, and 1% penicillin and streptomycin (Gibco, 
Grand Island, NY). Myoblast differentiation was induced as 
described before. 3637 On day 4 after induction of differen- 
tiation, cells were treated with pooled DMD serum (20 ul), 
or incubated with 5ng/ml TGF-p (Santa Cruz, Dallas, TX), 
10ng/ml FGF (Santa Cruz), or 50ng/ml TNF-a (gift from 
Dr.Yingying Le SIBS, CAS, China) for 12 hours before culture 
mediums (2 ml) were harvested. FGFR1 inhibitor PD1 73074 
was from Han-Xiang Chemical (final concentration 10 nmol/l). 

RNA isolation and qRT-PCR analysis. The miRNAs were 
extracted from serum and cell culture medium using miRcute 
miRNA isolation kit (Tiangen, Beijing, China). Small RNA was 
reverse-transcribed by using stem-loop RT which was per- 
formed according to the protocols described before 38 Real- 
time PCR was performed as described before 39 on an ABI 
7900 Real-Time PCR System (Applied Biosystems, Grand 
Island, NY). The primers for miR-1 33 detected both miR-1 33a 
and miR-1 33b. We chose miR-223 as an internal control for 
real-time PCR analysis of serum miRNA. Primers used in this 
study are listed in Supplementary Table S5. 

Myosin ATPase staining. Muscle fibers were counted, 
measured, and classified by ATPase reaction which was 
performed as described before with minor modification. 40 
Briefly, fiber types were determined by the differential stain- 
ing resulting from pre-incubation at either pH 10.5, 4.6, or 
4.2 for human fibers. The first method distinguishes type I 
(light staining) fibers from type II (dark staining) fibers, the 
second distinguishes type I (dark staining) from lla (light 
staining) and Mb or llc(medium staining) fibers, and the third 
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distinguishes lie (medium staining) from lla or Mb (bright 
staining)fibers. The three preincubation reagents were (i) 
20 mmol/l sodium barbital, 18 mmol/l CaCI 2 , pH 10.5, (ii) 50 
mmol/l sodium acetate, 30 mmol/l sodium barbital brought 
to pH on 4.6 with HCI, and the same as (ii) but adjusted 
to pH 4.2. The preincubation times were 15 minutes at pH 
10.5 and 5 minutes at the acid pH's. After preincubation, 
the sections were incubated for 30 mm in 20 mmol/l sodium 
barbital, pH 10.5, containing 9 mmol/l CaCI 2 and 2.7 mmol/l 
ATP; rinsed in three changes of 1% CaCI 2 (3 minutes each); 
immersed for 3 minutes in 2% CoCI 2 ; and rinsed in three 
changes of tap water. After staining for 1 minutes in 1% 
(NH 4 ) 2 S, the sections were washed with several changes of 
tap water, dehydrated with ethanol, cleared in xylene, and 
mounted in DPX. 

Statistical analysis. Data are presented as means ± standard 
error unless otherwise indicated. ANOVA and student's f-test 
were performed for statistical analysis as indicated. P values 
of < 0.05 were considered to be significant. ROC curves were 
plotted to evaluate the diagnostic effects of the profiles. The 
area under the ROC curve (AUC) was calculated for the mea- 
surement of discrimination accuracy. GraphPad software was 
used for all the statistical analysis. 
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Figure S1. Serum levels of myomiRs in dystrophic mouse 
model. 

Figure S2. The serum level of CK in BMD and DMD patients. 
Figure S3. The expression levels of miR-499, miR-208a, 
miR-208b, miR-206 and miR-21 in fast-twitch (tibialis ante- 
rior, TA) and slow-twitch (soleus, SOL) muscles of C57/B6 
mice. 

Figure S4. The impact of FGFR1 inhibitor PD1 73074 on 
myomiR levels in the medium of C2C12 myotubes. 
Figure S5. The effect of FGF, TNF-oc, or TGF-(3 on cell 
viability of C2C12 myotubes. 
Table S1. Patient information. 

Table S2. Correlation between age and serum miRNA/CK 
level of dystrophic children. 

Table S3. Correlation between age and muscle fiber type 
composition of dystrophic children. 
Table S4. Primer information. 
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